Deformed Gaussian Orthogonal Ensemble description of Small- World networks 



o 
o 

Oh: 

< 

in 



C/3 



I 

o 
o 



> 

o 
in 

oo 

o 



J. X. de Carvalho,i'B Sarika Jalan/'Q and M. S. Hussein2'[t| 

^ Max-Planck-Institut fiir Physik komplexer Systeme, 
Nothnitzer Strafe 38, D-01187 Dresden, Germany 
^Instituto de Ftsica, Universidade de Sao Paulo, 
CP. 66318, 05315-970 Sao Paulo, S.P., BraziE 

The study of spectral behavior of networks has gained enthusiasm over the last few years. In par- 
ticular, Random Matrix Theory (RMT) concepts have proven to be useful. In discussing transition 
from regular behavior to fully chaotic behavior it has been found that an extrapolation formula of 
the Brody type can be used. In the present paper we analyze the regular to chaotic behavior of 
Small World (SW) networks using an extension of the Gaussian Orthogonal Ensemble. This RMT 
ensemble, coined the Deformed Gaussian Orthogonal Ensemble (DGOE), supplies a natural foun- 
dation of the Brody formula. SW networks follow GOE statistics till certain range of eigenvalues 
correlations depending upon the strength of random connections. We show that for these regimes of 
SW networks where spectral correlations do not follow GOE beyond certain range, DGOE statistics 
models the correlations very well. The analysis performed in this paper proves the utility of the 
DGOE in network physics, as much as it has been useful in other physical systems. 

PACS numbers: 



INTRODUCTION 

Initiated by two seminal works [J, 0| , the last decade 
has witnessed a spurt in activities of network research 
Pi S S] ■ Regular and random networks are the two lim- 
iting cases of network topology. For the regular network, 
each node is connected in a fixed pattern to the same 
number of neighboring nodes; on the other hand, for the 
random network, each node is randomly joined with any 
other node. Real-world networks show the properties 
which are intermediate of the regular and the random 
one [E [Ij Si- ^'^^ example, many real- world networks 
from diverse field have very small diameter but have very 
high clustering, two characteristics shown respectively by 
random and regular networks. To model randomness and 
regularity. Watts and Strogatz proposed an algorithm to 
generate popularly known as Small- World (SW) network, 
which has the properties of small diameter and high clus- 
tering Moreover, this model network is very sparse, 
i.e. network with a very few number of edges, which is 
another property shown by real- world networks. 

The structure of networks is described by its associated 
adjacency matrix A. It is defined in the following way: 
Aij = 1 if i and j nodes are connected and zero otherwise. 
We consider only undirected networks. In this case, the 
adjacency matrix is symmetric and consequently has real 
eigenvalues. These eigenvalues give information about 
some basic topological properties of the underlying net- 
work Q. The fluctuations of these eigenvalues can be 
studied by Random Matrix Theory (RMT). 

There is a long history of applications of random ma- 
trix ensembles to model fluctuations of the spectra of di- 
verse systems ^]. Unfortunately analytical results exist 
only if some ideal conditions are fulfilled by the systems 
studied. On the other hand real physical systems usually 



depart from these conditions. In order to cover these sit- 
uations other ensembles have been introduced 01 • One 
such class of ensemble is the so-called deformed Gaussian 
orthogonal ensemble (DGOE) [1, [1, This ensem- 

ble has been proved to be particularly useful when one 
wants to study the breaking of a discrete symmetry in a 
many-body system such as the atomic nucleus. It is also 
useful for studying transition among classes of ensem- 
ble such as order-chaos (Poisson — > GOE) and symme- 
try violation (2G0E GOE) Recently Jalan and 
Bandyopadhyay show that spectra of various model net- 
works and real world networks follow universal random 
1^ . 13| , intermediate between Poisson 
Correlations among eigenvalues of 



matrix properties 



and GOE statistics. 
SW networks follow GOE statistics of RMT for certain 
range and after that they deviate from the GOE statistics 
We believe that the DGOE supplies a RMT basis for 
the Brody distribution and gives a more accurate de- 
scription of the GOE-Poisson transition than the Berry- 
Robnik [l^ model, which purports to justify the Brody 
formula from an RMT stand point. The Brody distribu- 
tion was used previously in SW statistics investigation 
l3, In the present paper we analyze the spectra of 
SW model networks using DGOE. Based on the results of 
reference IJ] we argue, and show through numerical sim- 
ulations that fluctuations of the spectra of the SW model 
follows the description of a transition Poisson-GOE. 



SMALL- WORLD NETWORKS 

Watts-Strogatz model of SW network is constructed by 
rewiring the edges of regular ring lattice with probabil- 
ity p. This rewiring procedure generates a network with 
some random connections, without altering the number 
of vertices or edges. For p = 0, structure of the regular 
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lattice or fc-nearest neighbor coupled network remains 
same; on the other hand, for p = 1, the regular lattice 
becomes random network. For the intermediate values 
of p, the graph is a SW network: highly clustered like a 
regular graph, yet with small characteristic path length 
like a random graph. This onset of SW property hap- 
pens for a very small value of parameter p. Characteris- 
tic path length is defined as the number of connections 
in the shortest path between two nodes, averaged over 
all pairs of nodes. For a network of size N and aver- 
age degree k, it scales as N/k if network is regular, and 
log{N)/log{k) if network is random. Clustering coeffi- 
cient (C) is defined as the ratio of connections between 
neighbors to the number of allowed links. For regular 
graphs C is very high (3/4), whereas for random graphs 
it scales as k/N. Small- world networks show intermedi- 
ate behavior between these two extremes, with average 
path length being as low as for the random graphs, and 
clustering coefficient as high as that of regular graphs. 
This intermediate statistical features of SW networks are 
reflected in their spectral fluctuations, and can be nicely 
described using the DGOE which provides a RMT basis 
for the deviation from the GOE behavior of the short 
range correlation aspect of the eigenvalues, exemplified 
through the spacing distribution, and the long range cor- 
relation measured by the A3 . In the following we supply 
a description of the GOE-Poisson transition within the 
DGOE. 



THE DEFORMED GAUSSIAN ORTHOGONAL 
ENSEMBLE (DGOE): TRANSITIONS AMONG 
UNIVERSALITY CLASSES IN RMT 



The joint probability distribution of elements of DGOE 
has the general form [9|] 



and variance 



P{H, a, (3) = Z^^ exp{-aTrH^ - f3TrHf ] 



(1) 



where Zn is a normalization factor and TrH is the trace 
of the matrix H. In order to describe two interpolating 
ensemble the matrix H must be chosen as the sum of two 
terms 



H = Hq + Hi, 



(2) 



where the matrices Hq and Hi define complementary 
subspaces of H. According to ([T]) for /3 ^ 00 the ele- 
ments of Hi vanish and H is projected onto the matrix 
Hq. Since in this work wc are concerned with the statis- 
tics intermediate between Poisson and GOE we will de- 
fine Hq as the Poissonian ensemble. It will be a diago- 
nal matrix with elements given by Hq.j = E^^iSij whose 
eigenvalues i?o.i are independent random variable with 
Gaussian distribution 
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(3) 



(4) 



The elements of the diagonal-less matrix Hi are also ran- 
dom independent variables with zero mean and variance 
given by 



/ tt2 \ _ 1 + ^ij _ \2 1 
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(5) 



where A = (1 + When /3 = (A = 1) the 

ensemble corresponds to the GOE. In the limit (3 ^ 00 
(A — 0), there will be only diagonal elements and the 
Poisson regime is attained. 
The average level density 



Xtt \N / Jq X 
and the cumulative level density 



(6) 



(7) 

were calculated by Bertuola et al. 17] , who observed that 
formula ([7]) provides a more accurate manner of unfolding 
the spectra than the usual polynomial unfolding used in 
[Tl ] . These formulas work very well in the regime close ei- 
ther to Poisson or GOE statistics. Intermediate between 
these statistics there is a transition regime characterized 
by a rapid change in statistics from almost Poisson to 
almost GOE. In this regime formulas ((6|) and ^ need 
corrections (see [isj|). 



SIMULATIONS AND RESULTS 

Numerical simulations of the SW networks are made 
by considering ensembles of 20 networks of size N — 2000 
and average degree k = 20. The adjacency matrix was 
diagonalized numerically and its first and last 300 eigen- 
values were discarded. Since an analytical expression 
for the average density is still lacking, the unfolding of 
eigenvalues was made by fitting the cumulative density 
or stair-case function 
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N{E)=Y,QiE-Ei 



(8) 



to Chebyshev polynomial using the linear least squares 
method. Ei are the eigenvalues of the SW network and 
O is the unit step function. 

For p = 0, the corresponding adjacency matrix would 
be a banded matrix with entries one in the band. As 
some connections are randomized with probability p, cor- 
responding adjacency matrix gets some non-zero entries 
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outside the band, at the expense of equal numbers of en- 
tries of one in the band. The mean value of the elements 
of these matrices is p and variance is Fig. (l)-(6) 

plot the adjacency matrix for different rewiring probabil- 
ities. Left sub-figure of (1) plots the adjacency matrix at 
the onset of SW transition {p ^ 0.002). For such a small 
value of p, very few connections are rewired and hence 
adjacency matrix is still almost banded with very few 
connections outside the band. Note that we take aver- 
age degree of network as k = 20, which leads to a sparse 
network (i.e. the number of connections is of the order 
of the number of nodes). Left sub-figure of Fig. (6) plots 
the adjacency matrix for p = 0.2, for this value of p, 20% 
of connections are rewired leading to the equal number 
of one outside the band. 

Results for the statistics intermediate between Poisson 
and GOE are obtained by diagonalization of an ensemble 
of random matrices. The mean value of the elements of 
these matrices were taken zero and the variance of the 
diagonal and off-diagonal elements given by ([4]) and ([5]). 
The unfolding of the spectra of the matrices is done using 
([7]) . In the simulations the values of a = 1 and the size of 
matrices N = 2000 are kept fixed. In order to simulate 
a transition Poisson-GOE, ensembles with 100 matrices 
and different values of A are considered. For each value 
of A we check between the density of eigenvalue given 
by ^ and the density of eigenvalues from the numerical 
calculation. If the agreement between the two is poor, 
the simulations are re-run using a corrected version of A, 
and a (called A and A in [l3|) These corrections 

are needed especially in the transitional regime alluded 
to following the discussion below Eq. (7). 

In the discussion of the deviation of the spacing distri- 
bution from that of Wigner, SW practitioners have used 
the Brody distribution, 15], which is given by. 



Pp{s) = As^exp{-Bs 



(9) 



where A and B are related to /3 through the normaliza- 
tion condition. Another distribution which also purports 
to describe the transition case was derived by [16| us- 
ing RMT and semiclassical considerations. The DGOE, 
which we use in this paper, supplies a natural RMT for 
the description of the GOE- Poisson and/or the Poisson- 
GOE transitions. 

In order to investigate the long-range behavior among 
the eigenvalues of SW model we use the Dyson-Mehta 
statistics A3. It is defined as 



A3(L;a) 



1 

— mm 

L Bi,B2 



a+L 



dE[x{E) - BiE - B2f (10) 



where Bi and B2 are obtained from a least-square fit. L 
is the average number of spacings in the integration in- 
terval, and x{E) the number of eigenvalues which are 
less than E (for DGOE x{E) is given hy ^ ). A3 
measures the least-square deviation of the function x{E) 



(the unfolded spectra) from a straight line in the interval 
[a, a -|- L]. In order to improve the statistics and avoid 
the introduction of correlations we choose successive in- 
tervals which overlap by L/2 [2^. According to RMT, for 
GOE the expected value for large values of L approaches 



A 
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3goe 



— [In L- 0.0687] 



(11) 



and for Poisson statistics it approaches L/15. 

In the following we present A3 results for SW networks 
for various p values, and corresponding DGOE. The near- 
est neighbor spacing distribution of SW networks, which 
probes for short range correlations of spectra, for the 
range of p < pc can be modeled by Brody parameter 
as described in [l^. After this values of p, which cor- 
responds to the SW transition as defined by Strogatz- 
Watts 1], the short range correlations of spectra still 
follows GOE statistics, but the long range correlations 
probed via A3 statistics follows GOE statistics only for 
certain range, and after that deviation from GOE statis- 
Which indicates possible breakdown 
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tics is seen 

of GOE theory for SW networks. And hence we turn to 
the random matrix theory of DGOE. Note that for p > Pc 
the range for which A3 follows GOE statistics depends 
upon the size and average degree of the network as well 
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FIG. 1: Left sub-figure plots adjacency matrix, and right sub- 
figure plots the spectral rigidity as a function of L for the 
adjacency matrix of SW model with p = 0.002 (circle) and 
the DGOE with A = 0.0065 (full line). A3 statistics for SW 
is plotted for 20 sets of random realizations of rewiring. 
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FIG. 2: Same as figure 1, for p = 0.005 and A = 0.0070. 
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FIG. 3: Same as figure 1, for p = 0.010 and A = 0.010. 
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FIG. 5: Same as figure 1, for p = 0.050 and A = 0.015. 
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FIG. 4: Same as figure 1, for p = 0.020 and A = 0.010. 
Also Poisson (dotted) and GOE(dashed) A3 are shown for the 
comparison. 



In Figs. (1) to (6), the spectral rigidity A3 are pre- 
sented for the different values of p. The values of p varies 
from pc = 0.002, corresponding to the onset of SW be- 
havior, to p = 1 which corresponds to a random graph. 
Each figure also depicts the A3 for DGOE describing a 
transition Poisson-GOE. The DGOE simulations are per- 
formed for matrices of size N=2000. Note that for each 
value of p it is possible to find a correspondent A such 
A3 for which DGOE fits A3 for SW model. The values 
of parameters p and A are listed in the table HI Using the 
criteria developed in we find that the critical value 
of A which separates the chaotic (random) from regular 
regime is Ac « 0.15. Therefore before this value of p 
the SW model is still in the regular regime, although the 
distribution of nearest neighbor spacing is totally com- 
patible with a GOE description. The multi-peaks in the 
density of eigenvalues for these values of p Hlj also 
supports this finding, because it indicates that the net- 
work still has large amount of regularity. In Fig. (3) to 
(6) the values for p are increased and the A3 comes into 
the chaotic (random) regime. As the value of p increases 
the spectra of SW becomes closer to the GOE prediction. 
In other words, the local regularity is gradually destroyed 
and the network becomes random. The DGOE descrip- 
tion which we are using to model SW to random network 
behavior, shows that for p > 0.05, behavior of A3 statis- 
tics can be modeled by a single value of A = 0.015. It 
suggests that under the framework of DGOE description. 
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FIG. 6; Same as figure 1, for p = 0.20 and A = 0.015. 

the network with p 0.05 has as much symmetry as for 
a complete random network {p 1). 

CONCLUSION AND DISCUSSION 

According to the RMT the Poisson statistic describes 
systems with localized states on certain bases and 
uncorrelated spectrum. On the other hand the GOE 
describes systems that become ergodic in the thermo- 
dynamic limit and have correlated spectra. For p = 
we have the ring graph which possesses TV symmetry 
(rotational symmetry). The numerical calculations of 
the spectra show several degenerate eigenvalues |l4| . 
There is no level repulsion and the spectra of ring graph 
should follow the Poisson statistics. However, as the 
value of the parameter p is increased gradually the 
rotational symmetry is destroyed and coupling among 
the eigenstates takes place. The spectra gradually 
suffer a transition from Poisson statistics to GOE. 
For p > Pc the spacing distribution, P(s), agrees with 
GOE description, however A3 statistics shows some 
part in the regular regime. This leads us to conclude 
that for p = Pc the local regular structure is destroyed 
and short-range correlation between eigenvalues is well 
described by GOE. However some residual local regular 
structure is still present and the long-range correlation 
among the eigenvalues measured by A3 is intermediate 
between Poisson and GOE. This residual regular struc- 
ture is merely connected to the symmetric nature of the 
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SW ring. This implies a symmetry constraint in the 
distribution and the existence of pseudo-periodic orbits. 
Such effects leads to a A3 which is a linear combination 
of a regular, term plus the GOE term Note 
that the P{s) is less sensitive to the finer details of the 
statistics than A3(L). The behavior of the SW level 
statistics ( both P{s) and A3(i)) in this regime can be 
completely modeled by DGOE which was constructed to 
deal with such situations ( Constrained GOE). Finally, 
for p > 0.05 the GOE description of P{s) and A3(L) is 
recovered. 

Before ending, we give a detailed assessment of the 
effect of the size of the random matrices on the results 
of the statistical analysis. We have extended our study 
above to sizes N = 500, 1000, besides N= 2000. For each 
case we have performed the simulations and the subse- 
quent DGOE analysis. Space limitation does not allow us 
to present our results in the form of figures but we have 
collected the relevant information in the table alluded to 
above. Tabled! 



p 


iV 


A 


A 




500 


0.0060 


0.0180 


0.002 


1000 


0.0034 


0.0116 




2000 


0.0065 


0.0845 




500 


0.0090 


0.0405 


0.005 


1000 


0.0050 


0.0250 




2000 


0.0070 


0.0980 




500 


0.0110 


0.0605 


0.010 


1000 


0.0065 


0.0422 




2000 


0.0100 


0.2000 




500 


0.0140 


0.0980 


0.020 


1000 


0.0085 


0.0722 




2000 


0.0100 


0.2000 




500 


0.0220 


0.070 


0.050 


1000 


0.0120 


0.144 




2000 


0.0150 


0.450 




500 


1.0 


500 


0.200 


1000 








2000 


0.0150 


0.45 




500 


1.0 


500 


1.000 


1000 


1.0 


1000 




2000 


1.0 


2000 



TABLE I: Results of the DGOE analysis of SW networks 
using different sizes of the random matrices. First column 
indicates the value of rewiring probability p, second column 
shows size, the third column is the DGOE transition param- 
eter, while the fourth column is the modified transition pa- 
rameter A — N\^. See text for details. 

The first column indicates the value of SW rewiring 
probability p which is allowed to vary from very small, 
0.002 to the allowed maximum of 1.00. In the second col- 
umn indicates the size of matrices. The last two columns 
indicate the deduced DGOE parameters A and A (see 
the discussion of the DGOE in the section following the 
Introduction). As a reminder, the parameter A, which 



takes the values inside the interval 0-1, measures the de- 
gree of deviation of the statistics from a pure GOE ( or 
pure Poisson) . The results shown in table H] clearly in- 
dicate that the SW network is a rigid GOE ensemble, 
regardless to the size for large values of p. The size does 
matter, however, for small values of p, where one sees a 
clear dependence of A on the size of the matrices used in 
the DGOE simulations. 

In conclusion, we have performed a statistical analy- 
sis of the SW networks within the DGOE. The analysis 
clearly demonstrates the usefulness of the DGOE statis- 
tics in supplying a solid basis of an RMT- based model to 
describe the chaos-order transitions in such networks. In 
general terms we conclude that there is a direct connec- 
tion between p and A, which points to a natural mapping 
of SW network onto the DGOE. Finally, for p = 1 when 
the system is totally random the GOE description is re- 
covered. 

From the random matrix point of view small-world 
networks studied here provide a very interesting system 
where depending upon the rewiring probability one can 
see that the short-range and the long-range correlations 
of the same ensemble of matrices belong to two different 
classes of random matrix models. From network point 
of view the analysis tells that on the one hand a small 
amount of random rewiring is enough to introduce short 
range correlations among eigenvalues suggesting spread- 
ing of randomness in the whole network, on the other 
hand DGOE statistics for long range correlations sug- 
gests the nature of symmetry in network. The future 
directions of this study is to understand the interplay 
of dynamical response which is based on the spec- 
tra of corresponding adjacency matrix and the symme- 
tries hidden in the network under DGOE framework. So 
far we have only concentrated on the small-world model 
network, providing a basis to the DGOE description of 
networks, future investigations would involve studies of 
real- world networks [25] . 



* Electronic address: 'josue@pks.mpg.de 

^ Electronic address: sarika@pks.mpg.de; Currently at : 
Department of Physics, National University of Singapore, 
2 Science Drive 3, Singapore 117542 

* Electronic address: hussein@if.usp.br 

S Supported in part by the CNPq and FAFESF (Brazil), 
and the Martin Gutzwiller Fellowship program at the 
Majc-Planck-Institut fur Physik komplexer Systeme- 
Dresden (2007/2008). 

[1] D. J. Watts and S. H. Strogatz, Nature 440, 393 (1998). 

[2] A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 

[3] S. H. Strogatz, Nature 410, 268 (2001). 

[4] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 
(2002) ; S. Boccalettia et al, Phys. Rep. 424, 175 (2006). 

[5] M. Doob in Handbook of Graph Theory, edited by J. L. 
Gross and J. Yellen (Chapman & Hall/CRC, 2004). 



6 



[6] T. Guhr, MuUer-Goeling e H. Weidenmiiller, Phys. Rep. 

299, 189 (1998). 
[7] F. J. Dyson, J. Math. Phys. 3, 1191 (1962). 
[8] M. S. Hussein and M. P. Pato, Phys. Rev. C 47, 2401 

(1993). 

[9] M. S. Hussein and M. P. Pato, Phys. Rev. Lett. 70, 1089 
(1993) 

[10] C. E. Carneiro, M. S. Hussein, and M. P. Pato, in Quan- 
tum Chaos, edited by H. Cerdeira, R. Ramaswamy, M. C. 
GutzwiUer, and G. Casati (World Scientific, Singapore, 
1991), p. 190. 

[11] J. X. de Carvalho, M. S. Hussein, M. P. Pato, and A. J. 

Sargeant, Phys. Rev. E 76, 066212 (2007). 
[12] J. N. Bandyopadhyay and S. Jalan, Phys. Rev. E 76, 

026109 ( 2007). 

[13] S. Jalan and J. N. Bandyopadhyay, Phys. Rev. E 76, 

046107 (2 007). 
[14] S. Jalan and J. N. Bandyopadhyay (submitted) 
[15] T. A. Brody, Lett. Nuovo Cimento 7, 482 (1973). 
[16] M. V. Berry and M. Robnik, J. Phys. A 17, 2413 (1984). 
[17] A. C. Bertuola, J. X. de Carvalho, M. S. Hussein, M. 

P. Pato, and A. J. Sargeant, Phys. Rev. E 71, 036117 

(2005). 

[18] A. J. Sargeant, M. S. Hussein, M. P. Pato, and M. Ueda, 
Phys. Rev. C 61, 011302 (2000). 

[19] There are a range of values of A where the formula ^ and 
(O, derived in [17| . must be revised. A way of doing this 
is to make the substitution: A'^ ^ A^{1 + (1 + A^/2)[l + 
^1 + AV4]} and a' a[l + + X''^l + Xy4]. 



[20] O. Bohigas, M. J. Giannoni, and C. Schmit, Phys. Rev. 

Lett. 52, 1 (1984). 
[21] Farkas, L J, Derenyi, Barabasi, A. L and Vicsek, T., 

Phys. Rev. E, 64, 026704 (2001). 
[22] M. S. Hussein, and M. P. Pato, Phys. Rev. Lett., 80, 

1003 (1998). 

[23] See, however, O. Bohigas and M. P. Pato, Phys. Lett. 
B595, 171 (2004). In this reference another effect was 
considered that leads to a A3(I/) which is a linear sum 
of a Poisson plus a GOE terms. If the statistical analy- 
sis is insensitive to, say, / percentage of the levels, then 
the long-range correlation measure is given by A3(L) = 
+ (1 - ffAsij^). For small values of /, which is 
the more likely case ( say, / = 0.1), the P{s) is hardly 
affected by this effect, but the A3{L) is significantly mod- 
ified, especially for large L, where the dependence is lin- 
ear. 

[24] M. Girvan and M. E. J. Newman, Proc. Natl. Acad. Sci. 
USA 99, 7821 (2002); L. M. Pecora and T. L. Carroll, 
Phys. Rev. Lett. 80, 2109 (1998); A. E. Motter, Y.-C. 
Lai and F. C. Hopensteadt, ibid. 91, 014010 (2003); R. 
E. Amritkar, S. Jalan and C. K. Hu, Phys. Rev. E 72, 
016212 (2005); C. Zhou, A. E. Motter and J. Kurths, 
Phys. Rev. Lett. 96, 034101 (2006); F. M. Atay, T. 
Biyikoglu, and J. Jost, IEEE Trans. Circuits and Sys- 
tems 1 53, 92, (2006). 

[25] J. X. de Carvalho, Sarika Jalan and M. S. Hussein (un- 
published). 



